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ABSTRACT 

Observations of metal lines in Lya absorption systems of small H I column density 
and their ubiquitous nature suggest that the intergalactic medium (IGM) was enriched 
to about Z ~ 0.01 Zq by a redshift z 3. We investigate the role of winds from 
small star-forming galaxies at high z in enriching the IGM. The existence of large 
numbers of small galaxies at high z follows naturally from hierarchical clustering 
theories (e.g. CDM). 

For analytical simplicity we assume that the galactic winds escape the galaxies 
at a single characteristic redshift Zin, and we model the galactic winds as spherical 
shock waves propagating through the IGM. We then calculate the probability distri- 
bution of the mctallicity of the IGM, as a function of time (for different values of Zin), 
adopting plausible galaxy mass functions (from Press-Schechter formalism), cooling 
physics, star-formation efficiencies, gas ejection dynamics, and nucleosynthesis yields. 
We compare this expected distribution with the observed distribution of metallicities 
in the Lya forest at z = 3, the metal poor stars in the halo of our Galaxy, and with 
other observational constraints on such a scenario. We find that galactic winds at 
high z could have enriched the IGM to a mean metallicity of Z ^ O.OIZ0 at z ^ 3, 
with a standard deviation of the same order, if Zin 5, and that this satisfies all the 
observational constraints. 
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1 INTRODUCTION 

Recent observations with the Keck 10 m telescope have 
shown that relatively low H I column density Lya absorption 
lines (A^Hi ~ 10^* cm~^) at z ~ 3 have corresponding metal 
lines. Cowie et al. (1995), Tytler et al. (1995) and Womble 
et al. (1996) reported that more than half of all Lya ab- 
sorption systems with A^hi ~ 10^'' cm~^ have CIV. These 
observations confirmed previous suggestions of enrichment 
in the Lya forest clouds (Meyer & York 1987; Tytler 1987). 

More recently, Songaila & Cowie (1996) have detected 
lines due to several other ions such as SIV. They have also 
claimed that the fraction of Lya lines with associated metal 
lines is larger than ^ 0.5 as reported earlier. They found 
C IV fines in 90% of clouds with Nr i > 1.6 x 10^^ cm'^ and 
75% of clouds with Nh i > 3.0 x lO" cm"^. They emphasize 
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that it would be difficult to estimate the exact fraction of 
Lya lines with associated metal lines with A^hi ^ 10^^ cm~^ 
as it would mean detecting C IV lines with Nc iv ~ 10^^ 
cm~^ which is difficult at present. It is also of interest that 
Songaila & Cowie (1996) found evidence of enhancement of 
the alpha process elements versus the Fe process elements, 
as is also the case in damped Lya systems (Lu et al. 1996), 
and in the metal poor stars of our Galactic halo. 

Such a large fraction of Lya lines with associated metal 
lines suggests that the intergalactic medium at 2: ~ 3 was 
enriched. The estimated metallicity in the Lya systems de- 
pends on the value of the ionization parameter assumed. 
Cowie et al. (1995), Womble et al. (1995) and Songaila 
& Cowie (1996) estimated the metallicity to be of order 
Z ~ 10~^ Zq for an ionization parameter in the range 
r ~ 10~^-^-10"^ ^ Ranch, Haehnelt & Steinmetz (1996), 
however, compared the result of their numerical simulation 
with the observed data and estimated a mean metallicity of 
Z ~ lO"^ '^ Zq. 
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The IGM could have been enriched in a number of 
ways, e.g., Population III stars, early- forming very massive 
(10^ to 10® Mq) objects (Carr, Bond & Arnett 1984) or 
galactic winds from galaxies at high redshifts. In this paper, 
wo investigate the possibility that supernova-driven galac- 
tic winds from small star-forming gala^xies at high redshift 
could have enriched the IGM. Such low-mass galaxies form 
naturally at high z in hierarchical clustering scenarios for 
galaxy formation, like CDM. If the gas in these galaxies can 
cool and form stars, then any metals produced when these 
stars go supernova can easily be ejected from the galaxy by 
supernova-driven winds, as the galaxies have shallow grav- 
itational potential wells (Dckcl & Silk 1986). This was in 
essence the scenario outlined by Couchmari & Rees (1986), 
who consider the subsequent influence of these galaxies on 
structure formation. Here we investigate this scenario in de- 
tail, incorporating up-to-date cooling physics and compare 
our models to recent observations, concentrating on issues 
most directly related to metal production. 

The amount of gas and the metallicity of the ejected 
gas from each galaxy depends on the mass of the galaxy and 
the parameters of star formation (e.g., Nath & Chiba 1995). 
The enriched gas would be carried in cooled shells which 
will propagate through the IGM, sweeping up matter. The 
shells (each with its own radius and metallicity, which have 
been fixed by the mass of the parent galaxy) will then collide 
and the enriched gas will be mixed. Given a long time, the 
gas will distribute more homogeneously and the resulting 
metallicity will depend on the mass function of galaxies, the 
redshifts when the galaxies erupted and the parameters of 
star formation (in addition to the cosmological parameters) . 

Several authors have commented on the production 
of metals as a possible by-product of astrophysical pro- 
cesses occurring during the early stage of galaxy format 
tion (e.g. Fukugita & Kawasaki 1994, Voit 1996). Specif- 
ically, if supcrnovae arc partly responsible for the gen- 
eration of the UV background (particularly at soft ener- 
gies; see Miralda-Escudc & Ostrikcr 1990, Madau & ShuU 
1996) and/or the reionization of the Universe (see Giroux & 
Shapiro 1996), then simultaneous cosmic metal production 
seems inevitable. Recently, Gncdin and Ostriker (1997) have 
presented the results of their numerical simulations of reion- 
ization and metal production at high redshift by Population 
III stars. MiraJdarEscude and Rees (1997) have discussed 
the production of metals by high redshift supernovae, but 
did not consider the details of the gas ejection dynamics 
and the mixing of the gas in the IGM. Even prior to these 
works. Silk, Wyse & Shields (1987) had contemplated on 
the enrichment of the IGM by bursting dwarf galaxies as 
a plausible cosmological process, continuing to the present 
day, along with the generation of the soft UV background. 

This paper is structured as follows. In Section 2 we out- 
line the chemical history of the Universe from the Big Bang 
to the present day, and place our model in this context and 
in the more general context of galaxy formation theories. In 
Section 3 we estimate the time evolution of the mean metal- 
licity of the IGM. We also estimate its variance, from the 
details of the collisions between the shells and the conse- 
quent mixing of the gas. For simplicity, wo assume a single 
epoch at which the galactic winds go off, assume the shocks 
to be spherical, and then model the evolution of the proba- 
bility distribution of the IGM metallicity. We then compare 



this distribution with the observed metallicities of the Lya 
clouds at z — 3, the metallicity distribution of the metal- 
poor stars in the halo of our Galaxy, and other observations 
in Section 4. 



2 THE CHEMICAL HISTORY OF THE 
UNIVERSE 

In this section we outline a plausible scenario that describes 

the chemical history of the Universe. We do this so as to 
place the model we present here in the wider framework 
of the galaxy formation problem and also to highlight the 
observational contexts in which our results are and are not 
valid. 

After the Big Bang, only hydrogen, helium, and trace 
amounts of lithium were present in the Universe; carbon and 
larger nuclei essentially did not exist. As the initial fluc- 
tuations in the dark matter spectrum grew, at some red- 
shift, the very high a peaks reached masses large enough 
that the baryonic component contained within these peaks 
may cool through molecular transitions (Blumenthal et al. 
1984, Couchman & Rees 1986, Tegmark et al. 1997). The 
relevant baryonic mass scales vary from 10^ to 10'' Mq for 
100 > z > 10 (see Figure 6 of Tegmark et al. 1997). The 
molecular hydrogen responsible for triggering the cooling 
mostly comes from H^-fH" interactions, where the H~ is 
made from H-fe~ reactions (a small number of free elec- 
trons remain after recombination; Couchman & Rees 1986). 
No metals exist yet to assist in the cooling (however, see 
Gnedin & Ostriker 1992 for a possible alternative). After 
collapse is initiated and the densities get high H-I-HH-H reac- 
tions presumably take over as the primary source of making 
H2. No dust grains are available so that absorption of H I 
onto dust, which is the dominant method of making H2 in 
star-forming regions in the Milky Way, is not important. If 
sufficient H2 is made, the cooled baryons may form stars. 
If this happens, then some of these stars might become su- 
pernovae and produce some metals and UV photons. But at 
very high z » 10, these high a objects are extremely rare, 
so that they do not produce significant amounts of metals 
(they may, however, reheat the IGM; Tegmark et al. 1997). 

As the age of the Universe increases, the much more 
common lower a peaks reach the mass scales required for 
the gas to cool, and baryonic collapse and possibly star for- 
mation becomes much more ubiquitous. It is this phase that 
we are concerned with in this paper. The mass scale here will 
be set mostly by atomic cooling, as molecular hydrogen will 
easily be photodissociated by the few UV photons produced 
by the few higher a peaks (Haiman, Rees, & Loeb 1996). The 
main uncertainty preventing a rigorous study of this is the 
stellar initial mass function (IMF) in these galaxies, where 
the collapsing gas is essentially unenriched. We know that 
it is probably more biassed to high masses than the local 
Salpeter IMF because otherwise at least 20% of halo stars 
would be pregalactic and made out of material that has only 
been processed once (i.e. during this phase; Miralda-Escudc 
& Rees, 1996). If this were true, then at least 20% of the halo 
stars would be heavy neutron-capture deficient (see Truran 
1981, also Section 4.2), which is unlikely (Ryan, Norris & 
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Beers 1996). In this work, we adopt a parameterization that 
takes this lack of knowledge about the IMF into account. 

Here wc adopt plausible galaxy mass functions (from 
Press-Schechter (1974) theory), cooling efficiencies, and nu- 
cleosynthesis yields to compute the amount of metals pro- 
duced. We then model the propagation of these metals by 
galactic winds through the IGM and compute the mean 
metallicity of the IGM and its spatial variance as a function 
of time. The mean value we get, about 0.01 Zq, compares 
well with the metallicities of the Lya clouds seen at z = 3, 
and wc conclude that the metals soon there could result from 
this phase of cosmic enrichment. In such a scenario, local star 
formation within the Lya forest clouds themselves does not 
contribute many of the metals. This would not be surpris- 
ing, as the systems are very small and diffuse (on the other 
hand, in the bigger, denser, and more metal-rich damped 
Lya systems, self-enrichment almost certainly happens). 

It is also probable that many of the lowest Z halo stars, 
in particular the Sr-doficiont stars, might have been made 
from material that was only processed during this phase (see 
Section 4.2). Most of the rest of the halo probably came 
from material that was enriched by subsequent generations 
of stars that formed either in smaller stellar systems that 
wore the Milky Way progenitors, or during the collapse of 
the halo (Eggen, Lynden-Bell & Sandage 1962). Another im- 
portant consequence of the enrichment process is the gen- 
eration of a large number of UV photons, which will raise 
the Jeans mass of the IGM, hence inhibiting further collapse 
of baryons into dark matter halos (e.g., Thoul & Weinberg 
1995). What this means is that the redshift range over which 
this phase happens may be quite small (in this work we ap- 
proximate it by a single characteristic value Zi„). Those UV 
photons may also in whole or part reionize the IGM (de- 
ponding on the efficiency with which they can escape the 
galaxies), therefore explaining the H I Gunn- Peterson effect. 

After this phase, star formation continues to happen in 
galaxies, although it may well have been supressed (par- 
ticularly in small galaxies; see Babul & Rees 1992, Efs- 
tathiou 1992) at z > 2 by the UV background produced 
by the galaxies described above, and by quasars at high red- 
shift. Large numbers of star-forming galaxies are observed 
in the field (see e.g. Broadhurst, Ellis & Shanks 1988; Cowie, 
Songaila & Hu 1991), but these are typically more massive 
than the galaxies that are relevant in the preceding two para- 
graphs (Cowie, Hu & Songaila 1995). Therefore any metals 
produced are unlikely to escape the galaxies because of the 
deep potential wells. Galaxies like the Milky Way have also 
experienced considerable star formation in the disk and have 
retained the metals produced; this is why the disk stars and 
the local interstellar medium (ISM) have metallicities close 
to solar, much higher than the 0.01 Z© produced in the phase 
above. Some of these galaxies may have since been tidally 
disrupted in interactions and the metals transferred from 
their ISM into the IGM. Because these galaxies might have 
had fairly metal-rich ISMs the IGM metallicity at z = 
might be somewhat higher than 0.01 Z©; strong observa- 
tional constraints on the z = IGM metallicity do not ex- 
ist at present. However, the Milky Way formed out of a gas 
cloud that was probably close to 0.01 Z© in mean metallicity 
(this is inferred from the peak of the halo metallicity distri- 
bution function of Ryan & Norris 1991), so that the IGM 



was not probably significantly enriched above this value by 
the time the Milky Way formed. 

The one case where the IGM metallicity at 2 = is well 
known, is of course in galaxy clusters, where the hot intra- 
clustcr medium (ICM) is metal-rich, with Z of 0.1 to 1 Zq. 
The measurement of the iron mass to optical luminosity ra- 
tio suggests that the cluster ellipticals provided these metals 
(Arnaud et al. 1992). The detection of Type II abundances 
with ASCA suggests that the metals come from Type II su- 
pernovae, which arc not observed to occur in present-day el- 
lipticals (Mushotzky et al. 1996, Loewenstein & Mushotzky 
1996). This suggests that these ellipticals might have under- 
gone a starbxirst phase early in their history (possibly an 
ultraluminous phase; see Sanders & Mirabel 1996), where 
substantial amounts of metals were made and which have 
since entered the ICM by cluster-related processes like ram- 
pressure stripping. 



3 GALACTIC WINDS AND THE IGM 

We will begin by explaining the different building blocks of 
our model — the galactic winds, their propagation in the 
IGM, and the interactions between the shells — and then 
compute the time evolution of the metallicity distribution. 

3.1 Galactic winds 

Galactic winds are believed to originate when the thermal 
energy of the interstellar gas in a galaxy exceeds its binding 
energy. The amount of gas and metals ejected by a galaxy 
of mass (including non-luminous matter) M depends on the 
parameters of star formation, viz., the IMF, star formation 
rate, and the efficiency of supernovae remnants in heating 
the ISM gas (Yoshii & Arimoto 1987). A variety of numerical 
works calculated the time for initiating the galactic wind and 
the metal enrichment as a function of these parameters (see, 
e.g., Matteucci & Gibson 1995 and references therein). We 
will use the results of Nath & Chiba (1995), which are based 
on simple assumptions, but are close to those obtained from 
more sophisticated approaches, for our purpose. 

Nath & Chiba (1995; hereafter NC95) estimated the 
fraction 7 of the galactic mass M ejected in a wind {Myji„d = 
7M) (their equations (10) and (14)), and the metallicity 
of the ejected gas, Z^i^d (their eqn (9)), using iron as the 
trace element. The total amount of metals ejected from each 
galaxy is then Z^i^dM-ujind- These values depend on the 
slope of the IMF {<j){m) oc m~^dm, x = 1.35 for Salpeter 
IMF) and the efficiency of the supernova remnants in heat- 
ing up the ISM gas (77 ^ 0.1). They incorporated the IMF 
slope, the upper and lower limits of masses of main sequence 
stars and the lower limit of the mass of stars which go super- 
nova, into a single parameter v^o, which was defined as the 
number of supernovae per 50 Mq of baryons. For a Salpeter 
IMF, X = 1.35, ~ 0.37 (for details see NC95). The im- 
portant concept here is that 7 and are single-valued 
functions of M given 1/50 and rj (typical examples of these 
functions are given in Figure 1 of NC95) . Consequently, the 
functions 7(M) and Zi„indiM), combined with the mass 
function of the galaxies, yield the total amount of metals 
produced. 
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The total energy in the explosion (meaning the cumu- 
lative effect of all the supernovae in a galaxy) is uncertain 
and we adopt the method of Tegmark, Silk, & Evrard (1993, 
hereafter TSE93) to estimate this. We need to know this 
number because it determines how well the winds and the 
metals they carry are propagated through the IGM. Since 
the galactic winds last for tb urn ~ 5 X lo"^ yr, which is much 
smaller than any other timescales in this calculation, we as- 
sume them to be ejected in an explosive manner. Making 
this assumption, the total energy in the galactic wind is 
then 0.02 x 0.007 x M^c^, where Mb is the baryonic mass in 
the galaxy. The 0.02 comes from the observations of Heck- 
man, Armus, & Milcy (1990) who show that 2% of the to- 
tal luminosity of high-z starburst galaxies goes into galactic 
winds. If one assumes the galactic wind to have a constant 
mechanical luminosity Lsn for the duration of tbum, then 
TSE93 show that Lsn ~ {Mb/MQ)LQ. In the calculations 
below, we assume Mb = 0.1 Af (we discuss the effect of the 
uncertainty in this relation on the final results in §4). 



3.2 Propagation of galactic winds in the IGM 

We will again follow TSE93 in calculating the propagation 
of the galactic winds in the IGM. The matter in the wind 
is assumed to be in a dense, cool shell (except for a small 
fraction of mass in the hot interior), that is sweeping more 
mass from the IGM. The density of particles in the shell 
can be estimated by equating the ram pressure {mpnaVX, 
where Ua is the ambient particle density and Vsh is the shock 
velocity) to the thermal pressure of the gas in the shell, 
which is assumed to be cold (at T ^ 10* K, the lowest 
temperature attained by gas in the absence of molecular 
cooling). This gives the particle density in the shell, = 
Tia 1201/3^ iQQ, where the shock velocity is written in the units 
of 100 km/s. Assuming the total mass inside the radius R 
of the shell to be contained within the shell, the thickness 
of the shell is 5R = (^/360)V;-J'ioo. 

Let the pressure inside the shell be p, which is driv- 
ing the shell outward with a force Ai^R^p. The sweeping 
of matter from the IGM provides a braking force equal to 
{R - HR)^ = {R- HR) x 3(1 - H), where H is the 
Hubble constant. There is another braking force due to the 
gravitation of mass in the shell. If the dark matter density 
(in units of the critical density) is Od, then the equation of 
motion of the shell is (eqn (1) of TSE93) given by 



R = 



_^^|__|(7?-//fl)2_(a,+o.5a,GM)(o.5H^i?).(i) 



Here, Qigm is the ratio of the density of the intergalactic 
gas to the critical density. This should be differentiated from 
the total baryon density in the universe, fib (expressed in the 
units of the critical density), eis some of the baryons would be 
in the collapsed halos and not in the intergalactic medium. 

The evolution of the total energy inside the shell {Et = 
IpV = 2npR^) is given by (TSE93) 



dEt 
dt 



—J— — Lsn — AixpR R — Lb 



(2) 



Here, the first term refers to the energy input due to super- 
novae, the second term is due to adiabatic loss, the third 
term describes energy loss due to free-free radiation and the 
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Figure 1. The uppermost panel shows the evolution of the radii 
of the shells from parent galaxies of mass M = 10* M0 (dotted) 
and M = lO" (solid), for 1 + = 5 and 7. The middle 

panel shows the corresponding shell velocities. The lowest pan- 
els shows the evolution of the porosity. The calculations assume 
sCDM model and Q.igm = 0.05. 



last term is the energy loss rate due to inverse Compton 
scattering off the microwave background photons. 

We solve the above two equations for the propagation of 
the shells in the IGM. For the initial conditions, we assume 
that the initial radius of the shell is the radius of the galaxy. 
We adopt the scaling of Saito (1979) for the typical size of 
a galaxy {R = 1.2 x 10^(M/10"AfQ)"-'''-' kpcj. The initial 
velocity is expected to be of the order of the thermal velocity 
of the ISM gas which is presumably heated to ~ 10® K 
(Babul & Rees 1992), i.e. Vsh,ioo ~ 1 initially. However, we 
found that the results arc almost independent of the initial 
velocity, as the kick due to Ls„ , which we assume to continue 
for 5 x lO'^ yr, basically determines the propagation of the 
shell. 

The upper two panels of Fig. 1 show the evolution of the 
radius and velocity of the shell ejected from galaxies with 
masses 10* and 10® M© at redshifts 1 -|- Zin = 5 and 7. 



3.3 Porosity 

If the comoving mass function of the parent objects which 
ejected their galactic winds at a certain redshift, Zin, is 
n{M,Zin) (which has the dimensions of Mpc~^ M^^) then 
the porosity por of the IGM at some later redshift can be 
defined as, 

--^por(M, z, Zi„) = ^R{M, z, Zi„)\{M, z,„) (-i±^) '.(3) 

dM 3 \l + Zin/ 

The porosity here simply denotes the filling factor of the 
spherical shocks (volume x number density of parent ob- 
jects) for por < 1. For larger values of por, the shells over- 
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lap one another, and in this case, the filhng factor of the 
overlapped regions is 1 — por~^ . The function n(M, Zin) can 
be estimated for a given model of structure formation from 
Press-Schechter theory. 

The Pross-Schochtcr mass function only refers to objects 
with a density contrast above a threshold value. To be spe- 
cific, it counts objects with linearly extrapolated overdensity 
5{1 + Zc) > 5c{l + Zc), where Zc is the collapse redshift. To 
form a galaxy by a given z the gas in the halo must, in 
addition to reaching some threshold density contrast, cool 
within a Hubble time tH{z). We take this effect into ac- 
count using the approach of Peacock & Heavens (1990). 
They showed that the effect of cooling can be taken into 
account by changing (1 + Zc) in the Press-Schechter mass 
function to (1 -|- z){l + M / Mcooi)^^^ , for galaxies to have 
formed by a given redshift z, where 

M,ool ~ 3.6 X 10" (i^) Mq ; (4) 
O.Uo 

here we use the nucleosynthesis value of ~ 0.05 h'^^ (ftso 
is the Hubble constant in units of 50 km s~^ Mpc~^) and 
we have assumed cooling in the absence of metals. We use 
this function for our calculations below. 

Formally, we should extend the above calculations to 
include the constraint that not only should the galaxy form, 
but it should also have had enough time for the supernovae 
to heat up the ISM and for the galactic wind to ensue. How- 
ever, we find that, for reasonable parameters of star forma- 
tion (i/^o = 0.35,77 = 0-1)) the ratio of the time scale for 
galactic winds {tgw, see NC95) and the cooling time of the 
gas tcooi is 



wind^wind 



M^ind + (47r/3)i?3pcfi/GM 



(6) 



tcooi V 109 Mo 7 V 6 ; 



(5) 



where we have assumed primordial cooling (Peacock and 
Heavens 1990), fij = 0.05 and hso = 1- The two time scales 
are of the same order for the range of masses we are inter- 
ested in. Therefore, the additional constraint {tcooi + tgw) < 
tn will not substantially change the above mass function, 
and so we neglect it here. 

As a template model of structure formation, we use the 
standard CDM (hereafter sCDM) model with /150 = 1 and 
n = 1. We use the normalization of as = 1.2 from the 
analysis of COBE data of four years (Bunn & White 1996). 
We use the analytical fit to the CDM power spectrum as 
given in White & Frcnk (1991) for our calculations. 

In the lower panel of Fig. 1, we plot the porosity as a 
function of redshift for M = 10* and 10^ M3 for 1 + z^n = 
5, 7 in the sCDM model, calculated by integrating eqn (3) 
with n(M') oc S{M' — M). The important thing to note in 
the figure is that shells from more massive galaxies reach 
the unit porosity limit later in time. 



3.4 Metallicity of the shells 

The metallicity of the shells will depend on (a) the amount 
of metals ejected by the parent galaxy {M^i„dZiaind', see 
above in §3.1), (b) the total amount of gas ejected {M^i^d), 
and (c) the amount of IGM gas swept up by the shell 
( (47r/3)i? Vcn/GM, assuming that it has not collided with 
any other shells so far; otherwise see below). The metallic- 
ity of the shell in this case is then given by 



If the enrichment of the IGM is absolutely homoge- 
neous, then the total metallicity of the IGM will be 



ZiGM, 



J M^i„dZwindn{M, Zin)dM 
J M^i„dn{M, Zin)dM + pcQiGM 



(7) 



Here the integration is over the range of galactic masses. 
We have used a lower limit of 10^ M© in the calculations. 
The value of ZiGM,hom will depend on the epoch of galactic 
winds Zin, the mass function of galaxies and the value of 
f2/GM, in addition to the parameters of star formation in the 
parent galaxies. For sCDM models, and for Q,igm = 0.05, we 
find that ^/cM.hom/^© = 0.033, 0.045, 0.057, for 1 + Zi^ ^ 
10, 7, 5 respectively. (These numbers do not change much for 
other values of O/gm; for Qigm = 0.01, for example, they 
are practically the same.) However, the state of the IGM is 
unlikely to be that with a constant metallicity, due to the 
inefficient mixing process. We next discuss the mixing due to 
collisions between shells. But before we do so, let us define 
another quantity, the homogeneous metallicity of the IGM 
due to galactic winds from galaxies of mass lower than a 
certain value M', 



ZlGM,hom{M' < M,Zin) 



M^indZwindn{M\ Zin)dM' 

J*^ M^indn{M', Zi„)dM' + pc^iiGM 



3.5 Interaction between shells 

When the shells meet, the result of each interaction will 
depend on the velocity of the shells involved. With the as- 
sumption of a single for all galactic winds, one can make 
a one-to-one correspondence between the velocity of the shell 
and the mass of the parent galaxy, and this makes it simpler 
to trace the evolution of the shells, as explained below. 

The interactions between the shells are very compli- 
cated in general. It will be useful here to divide them into 
two categories. The first type of interaction is that between 
shells of different sizes. For simplicity, let us assume that the 
interacting shells differ greatly in their sizes {Ri <^ R2, a,t 
some redshift z). In the approximation of a single epoch of 
galactic winds, as the curves in Fig. 1 show, the correspond- 
ing velocities are then also very different, with Vi <^V2. It 
is then also true (from the same curves) that the masses of 
the parent galaxies and therefore the metallicity of the shells 
are different: Zshell,i < Zshell,2- 

We assume that in this case, the shell with larger ve- 
locity and momentum (and therefore, with higher metallic- 
ity) will sweep the matter in the other shell with it. Let 
us look at the shells of size 7?2. Shells of much smaller size 
Ri <C R2 would have attained unit porosity much earlier 
than these shells. In other words, by the time the shells 
of size R2 have reached a reasonably large size, shells of 
much smaller sizes would have overlapped and distributed 
their metal enriched gas homogeneously in the IGM. For 
shells which have reached a size R2 at some redshift, the 
IGM gas already has a metallicity corresponding to shells of 
smaller sizes, originating from smaller galactic masses, that 
is, ZjGM,hom{M' < M2,Zin), where R2 is the size of the 
shell from a parent galaxy of mass M2. Although, strictly 
speaking, this is true only for R2 ^ Ri, for analytical sim- 
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Figure 2. The shell metallicity, Zghell.M, is shown as a function 
of the galactic mass M at different redshifts, for 1 + Zin = 10, 
sCDM model and njcM = 0.05. 

plicity we will assume it to hold also for R2 > Ri. We, 
therefore, modify our previous expression for the metallicity 
of the shells from galaxies with mass M (eqn 6) , and write 

ZsheU,M = 

MyjindZyjind + (47r/3)-R^n/GM pc ZiaM,hom{M' < M,Z, 

Myjind + {4w/3)R^pcniGM 

(For future reference, let us call this assumption Al.) This 
function Zsheii,M is a monotonically increasing function of 
M, a fact which will be useful to us below. 

We have seen (in §3.3) that shells from more massive 
galaxies reach the unit porosity limit later in time. Con- 
sider now the shells corresponding to mass (of parent galax- 
ies) M, . Once the porosity of these shells reach unity, the 
IGM gas, including the gas and metals from shells corre- 
sponding to M < Mt, has been more or less raked up by 
these shells. This means that ZsheH,M* does not increase any 
longer. Therefore, we freeze the metallicity of the shells once 
they reach unit porosity, and keep it constant during their 
evolution later. 

Fig. 2 shows the curves for Zsheii,M as a function of M 
at different redshifts, for zin = 10. In general, the shell 
metallicity decreases in time because of sweeping more and 
more of IGM gas, which must have lower Z. At later times 
(small z), the curves, therefore, flatten out. Note that, after 
reaching unit porosity, the shell metallicity does not change, 
giving rise to the envelope of the curves at the left hand side 
of Fig. 2. 

The second type of interaction takes place between 
shells of similar sizes {R\ ^ R2). This type of collisions in 
general gives rise to two reflected Shockwaves moving away 
from the plane of intersection, two contact discontinuities 
and another shock wave forming an annulus where the shells 
intersect. For simplicity, we neglect the annulus here. Fur- 
thermore, for geometrical simplicity, we assume that the re- 
gion bounded by the reflected Shockwaves can be reproduced 
by imagining the colliding shells passing through each other, 
since the reflected shocks have similar velocities as the origi- 
nal shells (see, e.g., Voinovich & Chernin 1995). The matter 
in this region is compressed and heated twice by the re- 



Figure 3. The two Umits of calculating the filling factor of gas 
of a given metallicity are shown schematicaUy. In Case I, the 
gas with the metallicity Z^i^^n m is confined within the shells 
originating from galactic mass M. In Case II, the gas occupies 
the overlapping regions. Note that the structure of shocks in the 
interacting region is very simplified here (see text). 

fleeted shocks to a temperature ~ 6 x 10^ K, for a typical 
shell velocity Vah ^ 50 km/s (see Fig. 1). The sound velocity 
of the gas at this temperature is ~ 30 km s~^. We assume 
that this gas fills up homogeneously the region bounded by 
the reflected shocks (or, in our oversimpliflcd picture, by the 
shells continued through each other). This is shown schemat- 
ically in the right panel of Fig. 3. Let us call this assumption 
A2. We note here that, this is strictly true only for shocks 
which do not have cold shell; in the case of cold shells, the 
gas forms a ring (Yoshika & Ikeuchi 1990). 

3.6 The probability distribution of metallicity 

Let us consider the galaxies of mass M. At some redshift 
z < Zin, the galactic wind shells from these galaxies have 
radii R. If the porosity for only those shells (as defined in 
§3.3) is less than unity (Case I in Figure 3), then these shells 
will have mostly collided with shells of smaller size, which are 
more numerous than shells of size greater than or equal to 
R. The shell will therefore have a metallicity Zsheii,M given 
by eqn (9) above (because of assumption Al). Since the gas 
is accumulated in the shell of thickness 5R, the filling factor 
of the gas of metallicity Zsheii,M is [{36R/R)por{M, z, %„)]. 
We could tentatively equate this as the probability of the 
IGM gas to have a metallicity Z = Z shell, m- However, be- 
cause of assumption Al, we note that the shells of size R 
will also be pushed around by shells of larger size, and there- 
fore, be mixed with higher metallicity gas (since Zsheii,M is 
an increasing function of M, or equivalently R; see Fig. 2). 
Strictly speaking, then, the above filling factor will equal 
the probability that the intergalactic gas has metallicity Z 
or higher. Let us call this quantity p{> Z). 

If the porosity is larger than unity (Case II in Figure 
3), then the filling factor of the gas with metallicity Z — 
Zsheii,M at redshift z is simply (1 — (l/por(M, z, Zi„))), from 
assumption A2. Again for the above reasons, we will equate 
this to p(> Z). 
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We can then consider galactic winds from galaxies of dif- 
ferent masses, calculate the shell mctallicities Z = Zsheii,M 
at each redshift, and calculate the function p{> Z) using 
the algorithm described above, depending on the value of 
the porosity for these shells. This function is essentially 
p(> Z) = p'{Z')dZ', where p'{Z') is the probability dis- 
tribution function of the IGM metallicity. 



4 RESULTS AND DISCUSSIONS 

We plot in the left panels of Fig 4a, 4b and 4c, the function 
p{> Z) at various redshifts for the cases 1 + Zin = 10, 7, 5. 
We assumed fi/GM = 0.05, vso = 0.37, 77 — 0.1, and sCDM 
model (with f2 = 1, hso = 1). The right panels show the 
the mean and the standard deviation as a function of the 
redshift for the three cases. The dashed line in the right 
panels show the values of ZjaM,hom, *-e., the metallicity of 
the IGM if the metals were distributed homogeneously. It is 
seen that the mean metallicity increases in time to approach 
this value. 

A note on the uncertainties involved in the calculations 
is in order here. We have done the calculations for other val- 
ues of the parameters Qigm and Lsn- Changing the value 
of fi/GM to 0.01, increases the shell mctallicities, and shifts 
the curves horizontally by at the most 0.3 dex. Decreasing 
Lsn by a factor of two, decreases the metallicities by ~ 0.15 
dex. 

4.1 Comparison with the metaUicity of Lya 

systems 

The reported metallicity of the Lya systems at « ~ 3 is 
Z/Zq ~ IQ-^-^-lO"^ (Rauch et al. 1996, Cowie et al. 1995). 
Fig 5 shows the mean metallicity and the standard deviation 
at 2; = 3 in our model, as functions of Zin- Curves are shown 
for fi/GM = 0.01 and 0.05. The figure shows that our model 
conforms to the observed values of the mean metallicity of 
the Lya absorption systems at a ~ 3 if ^ 5, considering 
both the observational and our theoretical uncertainties. 

It should be noted here that the determination of metal- 
licity in the Lya absorption systems depends on the assumed 
value of the ionization parameter F - the ratio of the den- 
sity of ionizing photons to that of the particles in the Lya 
systems which is uncertain. The resultant metallicity mea- 
sured is then uncertain to about a factor of two. It is hoped 
that future observations (including those of other elements 
than carbon) will give better values of the metallicity and 
its scatter, so that the above model can be tested against 
the data. In addition, abundance anomalies that arc seen 
in the Lya forest may yield important information about 
the masses of the stars that are producing the supernovae 
(Woosley & Weaver 1992), which in turn will allow better 
nucleosynthesis yields to be used in predicting the metallic- 
ities, in addition to information about F. 

4.2 Comparison with the metallicity of halo 
metal poor stars 

It is well known that there exist in the Milky Way halo a 
few stars with much lower Z than we would expect given the 




I^g.o VZs {1+z) 



Figure 4. (a) In the left panel, the function p(> Z) is shown 
at redshifts 2 = 7, 5, 3 for the case of 1 -|- Zin = 10. In the right 

panel, the mean (solid line) and the standard deviation (dotted 
line) of the IGM metallicity are shown as functions of redshift. 
The dashed curve shows the metaUicity had the gas from the 
galaxies distributed in the IGM homogeneously. The curves are 
for sCDM model and Qigm = 0.05. (b) The left panel shows 
p{> Z) at z = 5, 4, 3 for 1 -|- Zi„ = 7. The right panel shows the 
evolution of the mean (solid line) and the standard deviation of 
the IGM metallicity. (c) The left panel shows p(> Z) at z = 3 and 
2 for 1 + Zin = 5. The right panel shows the evolution of the mean 
(solid line) and the standard deviation of the IGM metallicity. 

metallicity of the initial gas cloud out of which the Galaxy 
formed. Six stars with Z < 0.001 Zq exist in the complete 

sample of Ryan & Norris 1991 (hereafter RN91 — see their 
Fig. 5d); the lowest metallicity star known is the carbon 
star G66-71, which has Z ~ IO'^^ '^Zq (Cass, Liebert & 
Wehrse 1988). The mean metallicity of halo stas is about 
Z ~ 10"^-*Zq (RN91). Such large inhomogcncitics might 
not have existed within the primordial gas cloud that was 
originally the Milky Way and are not required by this model 
(however they are not inconsistent with observation; RN91). 
The most natural interpretation is that these stars are pre- 
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Figure 5. The mean mctallicity of the IGM (solid lines) and the 
standard deviation (dotted lines) at z = 3 are shown as functions 
of Zi„, for iljGM = 0.01 and 0.05. The dashed line shows the 
reported metallicity of the Lya absorption lines at z = 3. 



galactic and dynamically are associated with the dark mat- 
ter, and not the rest of the stellar halo (however, binary 
mass transfer effects, like differential absorption of elements 
onto dust, might be important in explaining the low metal- 
licities of extremely low-Z carbon stars like G66— 71). We 
now investigate this possibility. 

A very approximate quantitative analysis is suggested. 
Although the probability distribution functions calculated 
above arc not Gaussian, let us discuss the case of a Gaussian 
distribution, for the sake of illustration. In our model, where 
the IGM metallicity has a mean and spatial variance of 0.01 
Zq, we would expect 18% of all pregalactic stars to have 
Z < 0.001 Zq. Here "pregalactic" refers to stars that formed 
out of material processed only through the phase described 
by this model and also to any stars that might represent the 
low-mass tail of the IMF describied in Section 3; stars that 
formed from gas that had been enriched further in smaller 
galaxies that were Milky Way progenitors are excluded. As 
6 stars in the RN91 sample satisfy this constraint, we would 
then estimate that 32 stars in their sample are pregalactic of 
which 16 would have Z < 0.01 Z©. This is 16% of the total 
number of stars (104) in the RN91 sample with Z < 0.01 
Zq, which in turn suggests 16% of the stars iu the halo with 
Z < 0.01 Zq are pregalactic. If one takes all the stars in the 
RN91 sample (372), then one finds 32/372 = 9% of all halo 
stars to be pregalactic. The rest of the stars formed during 
collapse of the halo, or in MW progenitors. 

So it is natural in our model that these stars should ex- 
ist. Our model also predicts the existence of a few very low Z 
stars like G66-71 (although in this case binary mass trans- 
fer is the more plausible explanation, particularly in light 
of the inverted C/0 abundance; Gass, Liebert & Wehrse 
1988). Surveys of the low Z halo stars, including both proper 
motion surveys like RN91 and objective prism surveys, are 
rapidly getting larger, and a metallicity probability distribu- 
tion function for the Galactic halo should soon be available 
with much better statistics at the low Z end than that in 
RN91. Then the calculation in the previous paragraph may 
be done rigorously. This will provide one of the strongest 



constraints on the model presented here, and more gener- 
ally on models that seek to explain the origin of the metals 
in the IGM. 

The pregalactic stars, as defined above, should have 
density profiles, like the dark matter, not r^"^ profiles, 
like most of the stellar halo (Rees 1997). Recent analysis 
of point-sources in the Hubble Deep Field (Elson, Santiago 
& Gilmore 1996) suggest that if the Milky Way has an ex- 
tended stellar halo, this extended halo must be very small 
compared to the rest of the halo. This suggests that the frac- 
tion of halo stars that are pregalactic is very low, and that 
our estimate earlier in this secion is too high. This prediction 
will not, however, be testable for some time because surveys 
of low Z stars would need to be carried out over regions well 
outside the solar neighborhood, where most stars are very 
faint. The best chance of making this measurement in the 
future is if we can identify a population of horizontal branch 
stars that are likely to be pregalactic based on metallicity 
and abundance ratio considerations. 

Another potentially powerful probe of the properties 
of pregalactic stars involves heavy neutron-capture element 
(e.g. strontium = Sr) abundances (see Truran 1981, also 
Magain 1989; Ryan, Norris, & Bessell 1991). In the sce- 
nario described by these authors, the first generation of stars 
(i.e. that described in this paper) only synthesizes iron peak 
nuclei from the unenriched gas. The second generation per- 
forms secondary neutron capture. Only third or later gener- 
ation stars would then have the heavy neutron-capture ele- 
ments. This would in turn mean that any primordial stars 
should be deficient in these elements. So would be any second 
generation stars, which are also among the ones we call "pre- 
galactic" here. The most detailed survey of Sr abundances 
comes from Ryan, Norris, & Beers (1996). Lower metallicity 
stars do tend to have less strontium. However, the evidence 
suggests that very few indeed, if any, are completely Sr defi- 
cient. Again, this means that the pregalactic fraction of halo 
stars that we estimated above is probably too high. It also 
means that the IMF in the galaxies in Section 3 is proba- 
bly biassed towards high masses (if it was Salpeter, then at 
least 20% of halo stars would have been pregalactic (specif- 
ically first-generation); Miralda-Escudc & Rees 1997). An 
important caveat is that accretion of material (in this case 
Sr) from the ISM onto the atmospheres of the halo stars is 
small (see Yoshii 1981), so that the Sr-deficient population 
we observe is the true primordial one. 

4.3 Constraints from the UV background 
radiation 

Gunn-Peterson tests show that the IGM was highly ion- 
ized by 2 ~ 5. Several models for the reionization of the 
IGM have been discussed in the literature, especially those 
investigating photoionization by various sources (quasars, 
star forming young galaxies, decaying neutrinos). Shapiro 
& Giroux (1987; see also Giroux & Shapiro 1996) claimed 
that ionizing photons from the observed quasar population 
were not enough for reionizing the IGM. Miralda-Escude & 
Ostriker (1990) showed that including the UV photons from 
the star forming galaxies at high redshift would make the 
spectrum of the UV background radiation softer. 

The UV background radiation at high redshift has been 
estimated using the proximity effect (Bajtlik, Duncan & Os- 
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triker 1988). It is known that the intensity of the UV back- 
ground at the Lyman limit is ~ ]^q-2i±o.5 ^^.g g-i 
Hz~^ sr~^ at 2 ~ 2-3, although the value at higher redshifts 
remains uncertain (see, e.g., Bechtold 1994; Cristiani et al. 
1996). Gunn- Peterson tests for both H 1 and He 11 atoms 
at high redshift give clues to the spectrum of the UV back- 
ground radiation. The recent observations of He H absorp- 
tion at X ~ 3 (see e.g. Songalia, Hu & Cowie 1995) have been 
used to further constrain the spectrum. Madau & Meiksin 
(1994) showed that the observed constraints allowed several 
possible sources, ranging from only observed quasars to a 
combination of quasars and star forming galaxies. 

Star formation responsible for the galactic winds at high 
redshift discussed above would necessarily produce UV pho- 
tons, and might reionize the IGM at high redshift. Silk, Wyse 
& Shields (1987) suggested that early galactic winds from 
dwarf galaxies would enrich the IGM and could also provide 
a soft UV background radiation. Miralda-Escudc & Ostriker 
(1990) noted that the intensity of the radiation is consis- 
tent with a cosmic metallicity of ^ 0.01 solar. Giroux & 
Shapiro (1996) considered the problem of metal production 
in galax;ies responsible for photoionization of the IGM in 
detail. They showed that the metal enrichment inside the 
galaxies could reach 20%-30% of solar metallicity by 2; ~ 3 
if the metals are not ejected from the parent galaxies. 

Madau & Shull (1996) considered the metal production 
rate in star forming galajcies at high redshift. They did not 
however consider physical processes that might distribute 
the metals outside the parent galaxies; their models dif- 
fer fundamentally from ours in that they assume that the 
metals observed in the Lya clouds at 2: = 3 were made 
inside those clouds. They considered continuous metal pro- 
duction in galaxies from z ^ 4 and determined the resulting 
UV background radiation at 2 ~ 3, after letting the UV 
photons traverse the cloudy IGM. For continuous star for- 
mation from 2 = 4, and for reaching a metallicity of 0.01 
solar in the parent galaxies by 2 = 3, the resulting UV back- 
ground radiation has an intensity (at the Lyman limit) of 
0.13 X 10~^^ erg cm~^ s~^ Hz~^ sr~^ This assumes that 
a fraction fesc = 0.25 of the UV photons escapes the par- 
ent galaxies. This value of the intensity is small compared 
to that estimated from proximity effect (see above). The 
important concept stressed by these authors is that fesc is 
very poorly known (see also Leitherer et al. 1995, and the 
discussion of that paper in Madau & Shull 1996, for obser- 
vational constraints). In addition, the IGM itself might have 
considerable opacity at 2 > 5 (no Gunn- Peterson constraints 
exist), and this might also reduce the contribution of the UV 
photons produced to the cosmic background. 

In our model, the galaxies emit their UV photons at 
higher redshifts than considered by Madau & Shull (1996). 
Given their calculations, the contribution of these photons 
to the UV background at 2 ~ 3 must be even less than pro- 
duced in their 3 < 2 < 4 continuous star-formation model. 
Given this, and the uncertainties above, it seems reasonable 
to state that the UV background we produce is likely to 
be significantly smaller than what is observed through the 
proximity effect. 



4.4 Can we see this process in progress? 

The brightest objects at Zin will be the supernovae that 
produce the metals here. They will have ~ 27 at 2 = 5 
(duration ~ 1 year) in the observer frame, which is too faint 
to see from the ground (Miralda-Escude & Rees, 1997). 

An interesting idea suggested by Miralda-Escudc & 
Rees (1997) is that a few of these high-2 supernovae will 
happen to fall behind rich clusters of galaxies and so have 
their total flux amplified by gravitational lensing. The cross- 
section that they estimate for magnifying a source by more 
than A = 100 is about lO"** arcminutes (this is the minimum 
magnification for such an object to be detected at present). 
This suggests that the probability of seeing such an object 
in any given rich cluster at a given time is 0.01%. If we 
surveyed 100 clusters to i<' = 22 every year, we would then 
have a 1% chance of finding one of these SN. As sensitive K- 
band detectors become available on large telescopes, lower A 
events become detectable and this project becomes feasible 
(optical detection is unlikely because observer-frame optical 
wavelengths probe rest-frame UV wavelengths, where super- 
novae are intrinsically faint). 

The discovery of such an object should be straightfor- 
ward to confirm spectroscopically, unless its 2 is so high 
(2 > 10) the strongest emission lines get redshifted too far 
into the infrared, where such faint-object spectroscopy is 
impossible. 

4.5 What is the fate of the small galaxies? 

In our model, the galaxies responsible for producing the met- 
als had their baryons blown out by the galactic winds (Dekel 
& Silk 1986). Therefore, if they survive to the present day, 
they must be almost dark, unless they have undergone re- 
cent infall of gas from the IGM and made more stars (this 
is the scenario described by Silk, Wyse, & Shields 1987). 

A few almost dark, low-mass galaxies arc known in the 
Local Group. The most extreme examples are Carina (Ma- 
teo et al. 1993), Draco, and Ursa Minor (Pryor & Kormendy 
1990) . Galaxies like Draco and Ursa Minor have stellar pop- 
ulations that are not self-gravitating and are > 99% dark 
matter by mass within their core-fitting radius. Their to- 
tal dark-matter fractions depend on the mass distribution 
outside this radius, and will probably be higher still. These 
galaxies are the natural remnants of the process described in 
this paper. The few stars that they do have probably result 
from later gas infall and subsequent star-formation. They do 
not have zero metallicity (Hodge 1989) so do not represent 
the low-mass tail of the IMF described earlier in this paper. 
Very few low-mass first-generation stars, if any, exist in these 
galaxies (a similar low-mass star deficient mass function has 
been suggested for the Galactic halo; Nakamura Kan-ya & 
Nishi 1997). 

But these galaxies do not exist in largo number. The 
galaxy luminosity function in the Local Group (van den 
Bergh 1992) is certainly not increasing rapidly at Mb ~ —7, 
the luminosity of Draco & Ursa Minor. 

So where have the galaxies that produce the metals 
gone? Many might have undergone considerable subsequent 
gas infall (as argued in Silk, Wyse & Shields 1987) and so be 
luminous galaxies today. They may even have merged with 
larger galaxies, or even have been tidally disrupted. Draco, 
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Ursa Minor, and Carina, if they really are remnants of the 
present model, would then be anomalous in that they had 
not suffered any of these fates 

An alternative is that more recent gas infall has not 
happened in most remnants of the present process, and that 
there exist many undetected dark galaxies today. Draco and 
Ursa Minor would then be anomalously bright examples 
of these objects, where the little amount of gas they picked 
up allowed them to form a few stars, which in turn allowed 
them to be detected. 



4.6 The origin of the X-ray emitting gas in 
clusters 

Presumably the X-ray emitting gas that is now seen in clus- 
ters has experienced the process described in this paper (al- 
though the galaxies here are not the source of most of the 
ICM metals; see Section 2). In two previous papers (Tren- 
tham 1994, NC95), we argued that dwarf galaxies might 
contribute much of the ICM. In the context of the present 
model, much of the IGM might well have been processed in 
small galaxies, but these are very different systems from the 
present-epoch cluster dwarf galaxy population discussed in 
those papers. 



4.7 C Concluding remarks 



To summarize, we have found that galactic winds at high 
z could have enriched the IGM to a mean metallicity of 
Z ^ O.OIZ0 at 2: ~ 3, with a standard deviation of the same 
order, if Zi„ ^ 5, and that this satisfies all the observational 
constraints. 



Wl lile our results show that these kinds of processes can 



easily account for the observed IGM metallicity, we needed 
to make a number of simplifying assumptions and adopt con- 
venient but plausible values for our model parameters. This 
work has been presented mostly in the spirit of a plausibil- 
ity argument, and we have not attempted a detailed study 
of parameter space. Numerical investigations will hopefully 
provide more insight into this process and will allow an as- 
sessment of the permissible ranges for the input parameters 



in models like the present one. 'i'hese will provide important 
constraints on galaxy formation tneories. 

Additional observations will also provide stronger con- 
straints. Specifically, measurements which should be partic- 
ularly valuable are (i) measurement of abundances ratios in 
the Lyof forest, (ii) the metallicity distribution function of 
the Galactic halo for a large sample of stars, and (iii) heavy 
neutron-capture elemental abundances for an equally large 
sample. 

B.N. thanks Drs. M. Chiba, P. Biermann and N.T. 
thanks Drs. L. Cowie, J. Norris, and M. Rees for stimu- 
lating discussions. The authors also thank the anonymous 
referee for his comments. 
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